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Abstract. We consider a particular class of n-dimensional homogeneous dif- 
fusions all of which have an identity diffusion matrix and a drift function 
that is piecewise constant and scale invariant. Abstract stochastic calculus 
immediately gives us general results about existence and uniqueness in law 
and invariant probability distributions when they exist. These invariant dis- 
tributions are probability measures on the n-dimensional space and can be 
extremely resistant to a more detailed understanding. To have a better anal- 
ysis, we construct a polyhedra such that the inward normal at its surface is 
given by the drift function and show that the finer structures of the invariant 
probability measure is intertwined with the geometry of the polyhedra. We 
show that several natural interacting Brownian particle models can thus be 
analyzed by studying the combinatorial fan generated by the drift function, 
particularly when these are simplicial. This is the case when the polyhedra is 
a polytope that is invariant under a Coxeter group action, which leads to an 
explicit description of the invariant measures in terms of iid Exponential ran- 
dom variables. Another class of examples is furnished by interactions indexed 
by weighted graphs all of which generate simplicial polytopes with n! faces. 
We show that the proportion of volume contained in each component simplex 
corresponds to a probability distribution on the group of permutations, some 
of which have surprising connections with the classical urn models. 

1. Introduction 

Consider the following two examples of stochastic processes. The first one, called 
the bang-bang process is classical and is particularly important in stochastic control 
theory. See the articles by Karatzas and Shreve [12] , Shreve [15] , or the chapter by 
Warnecke [16] . This is a one dimensional diffusion which solves the following SDE 
with a single real parameter a: 

dX t = -a sign(X t )dt + dW t . 

It can be easily shown that the process is recurrent if and only if a is positive. In 
that case, the process has a unique reversible stationary distribution under which 
\Xt\ is distributed as Exponential with rate 2a. 

The other example is the class of Brownian motions with rank-based interactions. 
This is a family of n one dimensional diffusions which is parametrized by a single 
vector S in W 1 . These diffusions have an identity diffusion matrix and a drift that 
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depends on the order in which the coordinates can be arranged in increasing values. 
If we think of each diffusion as recording the position of a particle moving on a line, 
then at any instant of time the particle with the ith smallest position gets an 
instantaneous drift Si. The formal SDE for the diffusion can be described by 

n 

(1) dX t (i) = Y,$A{Xt(i) = X t \j] }dt + dW t (i), i = l,2,...,n, 

where X t [l] < X t [2] < ... < X t [n] are the coordinates arranged in increasing 
order. The Wt(i)'s are asumed to be independent Brownian motions for some 
suitable underlying filtration. 

The rank-based interacting Brownian motions or closely related models have ap- 
peared in several veins of the literature. Extensive reviews can be found in the 
articles by Pal and Pitman [T3] and Chatterjee and Pal [3]. Some of the recent 
work include the articles by Jourdain and Malrieu [llj and Banner, Fernholz, and 
Karatzas [2]. Also see the related discrete time models by Ruzmaikina and Aizcn- 
man [14], and L-P. Arguin pQ. 

The recurrence, transience, and other long term properties of the rank-based 
Brownian motions can be precisely determined. The following result is from [13] . 

Theorem 1 (Theorem 4 in [13)). For 1 < k < n let 

k i n 

(2) a k :=J2(6i-S), S:=-J2$i- 

i=l i=l 

For each fixed initial distribution of the n particle system with drifts Si , the collection 
of laws of Xt[n] — Xt[l] for t > is tight if and only if 

(3) a k > for alll < k< N- 1, 

in which case the following result holds: 

The distribution of the spacings system {X t [j + 1] — A t [j], 1 < j < n — 1) at 
time t converges in total variation norm as t tends to infinity to a unique stationary 
distribution for the spacings system, which is that of independent Exponential vari- 
ables Yk with rates 2ctk, 1 < k < n— 1. Moreover, the spacings system is reversible 
at equilibrium. 

The independence of the spacings under the invariant distribution is somewhat 
puzzling since, due to the interaction, there is no independence between the spac- 
ing processes. The proof, which does not shed light on this phenomenon, invokes 
Williams's results on reflected Brownian motions [17] . 

There are some similarities between the two examples. Both involve drift func- 
tions that are piecewise constant; in fact the drift is a single vector under the action 
of a group (sign flips for the former and permutations for the latter). Moreover, 
in both cases the invariant distribution involves independent Exponentials which 
provide a friendly description of an otherwise abstract probability measure. 

We provide alternate proofs of these results as an application of the following 
general theory. Notice that the drift function in both the examples is the negative of 
the gradient (in the sense of distributions) of a positively homogenous (homogenous 
of degree one) function in M. n . For the bang-bang process, this function is a\x\, 
while for the rank-based processes it is given by — Y^i=i^ iX l^\- continuous 
positively homogenous function which is nonnegative for all x can be viewed as the 
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Minkowski (or, the gauge) function of a certain body containing the origin. Since 
we consider piecewise constant drifts, these bodies are n dimensional polyhedras. 
The corresponding stochastic process is shown to be recurrent when the polyhedra 
is a bounded polytopcs, which in turn happens when k is strictly positive for all 
i^O. When this is the case, the invariant distribution for the diffusion can be 
obtained from the uniform distribution on the polytopc. 

Let us now focus on the uniform distribution on an n-dimcnsional polytopes. If 
the polytope is a unit simplex {x : x% > 0, y]^ Xj < 1}, the uniform distribution 
can be effectively generated by dividing iid Exponential random variables by their 
total sum. Now suppose a polytope can be triangulated, i.e., decomposed as a 
union of K simplices with disjoint interiors. It is a well-known theorem that such 
a triangulation is possible for every polytope. Every component simplex, say Si, 
in the triangulation is a linear map Ai of the unit simplex. Thus, an algorithm 
to generate a point uniformly from the polytope would be to generate a point X 
uniformly from the unit simplex, choose / between {1,2,..., K} with probability 

(4) P(I = i)= J° l(Si) , i = l,2,...,K, 

Ef=iVol(^-) 

and let Y = A]X. Such an Y is clearly uniformly chosen from the polytope. 

A particularly explicit triangulation exists when the polytope is simplicial. That 
is, each of its extremal face is an (n — 1) dimensional simplex. One can then 
simply connect the origin (which is in the interior) with each of these faces to 
generate a nice triangulation. How does one check if a polytope is simplicial ? 
We demonstrate a simple condition when the symmetry group of the polytope is 
Coxeter, i.e., generated purely by reflections. Coxeter groups find applications 
in several ares of mathematics. The finite Coxeter groups include symmetries of 
regular polytopes and the Weyl groups of simple Lie algebras. They are usually 
defined formally as a set of generators and relations among them. However, we 
consider them in their original geometric form as treated by H. S. M. Coxeter in 
the classics [3] and [5]. The definitions and properties of irreducible group actions 
and Coxeter groups have been described in Subsection 13.11 This is a particularly 
nice case, when not only the polytope is simplicial, but it is also regular, and hence 
the random variable I is uniformly distributed among {1,2,..., K}. To connect the 
dots with independent Exponentials, we simply need to describe the maps A^s. In 
fact, due to regularity, determining A\ is enough, since the other maps are merely 
orbits under the group action. This is exactly the case for the bang-bang or the 
rank-based processes. The regularity corresponds to exchangeability among the 
particles, i.e., the condition that if the initial distribution is exchangeable among 
the coordinates, then so is the distribution at every other point of time. 

The second case we consider is not regular and does not involve any groups. The 
interaction is parametrized by all graphs with n vertices and possible edge- weights. 
In this case, the maps ^4i's are simple and explicit. However, the probabilities in (j4|) 
are not. In fact, these probabilities correspond to probabilities of various orderings 
of particles, increasingly arranged, under the invariant distribution. Hence, these 
polytopes induce probability distributions on permutations of n labels. We take up 
a few examples and show surprising connections with existing probability models 
on permutations. 

In the regular case, one of the results we prove is the following. 
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Proposition 2. Consider the SDE 

(5) dX t = b(X t )dt + dW t , where W t = (W t (l),W t (2),...,W t (n)) 

is an n-dimensional Brownian motion. Assume that the drift junction b : W l — » M" 
is piecewise constant and satisfies 

(6) b(ax) = b(x), V a > 0, and, b(Ax) = Ab(x), V A G G, 

where G is a finite irreducible group of orthogonal matrices. 
Then the following conclusions hold. 

(1) Le£ fc(ai) = — (x,b(x)). A sufficient condition for X to be recurrent is that, 
for some non-zero vector X G M n , we have 

k(x) = max (AX, x) . 
AeG 

In that case X has a unique, reversible invariant distribution [i on M" . The 
marginal law of X t converges in total variation to fx as t tends to infinity. 

(2) // G is Coxeter, there is a set of n many linearly independent vectors 
{Vii Vz-i • • • i Vn} o-nd n many constants {ct\, &i, . . . , a n } such that under /i, 
the random variables 

Yi = cti(Ar]i,x) , if k(x) = (AX,x), i=l,2,...,n, 

are iid exponential random variables with rate two. 

(3) Additionally, if the stabilizer subgroup of X in G is trivial, then the vectors 
{jlXi f]2t ■ ■ ■ j Vn} cire determined as the generators of the conic hull of the 
set of vectors {X — AX, A G G}. That is to say, every vector in the set 
{X — AX, A G G} can be represented as a linear combination of the subset 
{r)i, r)2, . . . , T] n } with nonnegative coefficients. 

The constants ct\,...,a n are the unique positive coefficients of r\\, . . . , r\ n 
in the expansion 

n 

X = aiT h- 
i=l 

Note that, we really do not need to know the details of the group structure to 
apply the previous result, except for the information that G is Coxeter. As we 
show in the examples, necessary and sufficient conditions can be obtained if wc 
have a better knowledge of the group structure. Finally, let us mention that a list 
of Coxeter groups up to isomorphisms is available and can be found in any standard 
textbook, say [10] . 

In Subsection 13.21 we describe several families of interacting diffusions that can 
be analyzed by the previous theorem. They all appear as solutions to stochastic dif- 
ferential equations of the type ([5]) with a piecewise constant drift function satisfying 
conditions ([5]), but involving different families of orthogonal groups. 

When the group is the group of permutation matrices, we get back rank based 
interactions. Using Proposition [21 wc provide an alternative proof (not involving 
reflected Brownian motions) of the Pal-Pitman result. 

The second class of examples are called sign-rank based interactions. Here the 
drift vector not only changes when the coordinate values get permuted, but also, 
when when the signs of coordinates change. The relevant group is the one which 
generated by all the permutation matrices and all the diagonal matrices whose 
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diagonal elements are either plus or minus one. In one dimension, this boils down 
to the simple Bang-bang process. 

The third class of examples are similar to sign-rank based processes, but with 
more constraints. Here, too, the drift vector changes when we permute coordinates. 
It also changes when we change signs of coordinates, but only when done in pairs. 
The group behind the curtain is generated by permutation matrices and diagonal 
matrices whose diagonal elements are ±1 with the additional constraint that only 
even number of — l's are allowed. 

Readers acquainted with the theory of Coxeter groups will recognize that the 
previous three examples correspond to the three well-known families of Coxeter 
groups, denoted by A n , B n , and T> n for each n G N. In each case Proposition [2] 
allows us to formulate a simple sufficient condition for checking the existence of a 
unique invariant probability distribution and provides a complete description of the 
distribution in terms of independent Exponentials. 

In the case of interaction through graphs wc consider the following class of inter- 
acting diffusions. Let Q be a graph on n vertices where the vetices are labeled by 
{1,2, ... ,n}. The edge between i and j have an associated edge weight fy, which 
is zero if there is no edge between the two vertices. Consider the SDE on R™ given 

by 

n 

dX t {i) = Assign (X t {j) - X t (i)) dt + dW t (i), i = l,2,...,n, 
j=i 

where, W is again an n-dimensional Brownian motion. 

When all the edge weights are nonncgativc, the model can be described by saying 
that the Brownian motions, which are indexed by the vertices of the graphs, get 
attracted towards one another. The constants (3ij measure the strength of their 
attraction. 

Unless the graph is the complete graph with constant edge-weights, the interac- 
tion is not regular. However, if we define 

1 " 

i=i 

the centered vector (X(l) — X, X(2) — X, . . ., X(n) — X) have an invariant distri- 
bution whenever the graph is connected and the edge-weights are nonnegativc. 

For this class of interactions our main focus of investigation is the law of random 
permutation that takes indices of coordinates to their ranks under the invariant 
distribution. This is not uniform by virtue of not being regular. As an interesting 
example, we consider the case when we assume that each particle has a mass rrii, 
i = 1,2,..., n, and that fyj = mirrij in the sense that the strength of the mutual 
attraction is proportional to the product of their masses. Due to this gravitational 
intuition, under the invariant distribution, wc should expect heavier particles to 
stay at the middle of the pile, while the lighter ones should be at the edge (the Sun 
being at the center, and Pluto at the far end). 

In general, this is very difficult to prove. However, in one particular case, this 
becomes apparent. For any a > 0, consider n particles with the interaction de- 
scribed in the previous paragraph, where the mass of the first particle is a and the 
rest of the masses are 1. Let <r(l) denote the rank of the first particle under the 
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invariant distribution of the centered vector (X(l) — X, X(2) — X, ... , X(n) — X). 
A surprising connection with Polya's urn scheme emerges. We prove the following. 

Proposition 3. Consider a Polya's urn scheme which initially has 2a red balls 
and 2a black balls. At every step one picks ups a ball at random, returns the ball to 
the urn and adds an extra ball of the same color. Then, the distribution of cr(l) — 1 
is the same as the number of red balls picked when we run the urn scheme described 
above for n — 1 steps. 

In particular, the sequence of random variables o~(l)/n converges weakly to the 
Beta(2a, 2a) distribution as n tends to infinity. 

1.1. Outline of the paper. In the next section we describe the set-up of the paper 
and prove general results about recurrence of interacting diffusions and their invari- 
ant distributions when they exist. Section [3] describes the combinatorics involved in 
the invariant distributions. Subsection 13.11 describes the connection with Coxeter 
groups followed by several examples in Subsection 13.21 The following subsection 
13.31 proves results about interactions parametrized by graphs. 

2. Diffusions with piecewise constant drift 

Consider a sequence of n-dimcnsional cones C\ , C2 , ■ ■ ■ , C r whose interiors are 
disjoint and the closure of their union is the entire space R™. Ignoring their mutual 
intersections (a set of measure zero) they can be thought of as a partition of R". 
Let b : R™ — * R n be a function that is constant over each Cj. In this section we 
prove some general results about the class of diffusions which satisfy the following 
stochastic differential equation: 

(7) dX t = b{X t )dt + dW t , 

where W t = (W t (l), W t {2), . . . , W t (n)) is an n-dimensional Brownian motion. 

The existence and uniqueness in law is immediate by an application of Girsanov's 
theorem. Define the function k(x) = — (x,b(x)). Since b is constant over cones, 
it follows that k is a positively homogeneous function. For example, k could be 
the Minkowski functional (i.e., the gauge function) of a convex body containing the 
origin. If k is continuous, by virtue of being piecewise linear, it follows easily that 
b is the negative of the gradient of k in the sense of distributions. In that case the 
SDE in ([7]) is an example of the Langevin equation. The following proposition is 
well-known about Langevin SDE. 

Proposition 4. Consider the stochastic differential equation ([7]). Let k : R" — > R 

be a continuous function such that b represents — Vfc in the sense of distribution. 
Assume that exp(— 2k(x)) is integrable. Then the probability distribution given by 
the un-normalized density exp(— 2k(x))dx provides a reversible, invariant probabil- 
ity distribution /i for the process in ([7]). 

Proof. Let \x be the measure on (R n ,£(R n ) defined by fi(dx) = e~ 2k ^dx. Consider 
the Sobolcv space, if 1 ' 2 , of all measurable functions / : R" — > R, such that / and 
all of its partial derivatives df/dxi, i = 1,2, ... ,n (in the sense of distributions) 
arc in L 2 (/i). Then wc can define the following symmetric bilinear form, on the 
domain H 1,2 , given by 

£(/,<?)= / (V/,V. 9 )e^W 
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Since e~ 2k ^ is never zero, it follows that H 1 ' 2 is a Hilbert space. Thus, £ is closed, 
since it is defined everywhere on the Hilbert space H 1 ' 2 . It is also known to be £ 
is Markovian (sec, e.g., |i7i|, example 1.2.1.). Thus, it is clear that this is a Dirichlct 
form in L 2 (R",/i). 

By Theorem 1.3.1 in [7], we claim the existence of a unique non-positive definite 
self-adjoint operator C! on H 1 ' 2 such that 

£(/,.g) = (v /z Z 7 /,v / ^C 7 . 9 ), t , V/, ff e 

Here (•) refers to the usual inner product in L 2 (/i). Or, in other words, (Corollary 
1.3.1 of [7]) there is a unique self-adjoint operator £' on a domain T>(C) C H 1 - 2 
such that 

(8) £(f,g) = (-jC'f,g}^ V / G £>(£'), V g G H 1 - 2 . 

We now show that ||5J) is satisfied by a multiple of the generator of the Markov 
process in ([7]). The generator, £, is given by 

Cf = (b,Vf) + ±Af. 

By our assumption b = — Vfc takes finitely many values. Thus, we can define C on 
the domain 



H 2 > 2 = {fe L 2 (/i) -L g L 2 (//),and —4 € L 2 (xt) Vi = 1,2. 



<9 2 f 

^-eL 2 (,),and^ 

It is clear that the domain of C above is a subset of H 1 ' 2 . 

We claim that 2£ satisfies (jHJ. In that direction, consider any / g if 2,2 and any 
g g Jf 1,2 , we have 

C{f)ge-^dx = ((b,Vf) + ±A/) ge'^da 

(b,Vf)ge- 2k ^dx+ [ ^Afge-^dx 
/ (Vk,Vf)ge- 2k ^dx-l [ (Vf,V(ge- 2k W))dx 
= -\J WMe-Wdx = -\£{f,g). 

Note that the boundary terms are zero in the integration by parts above since both 
df / dxi and g are in L 2 (/x) , and thus 



(9) 



^L ge -2k(x) 

dxi 



= 0, Vi = l,2,...,n. 



We can rewrite (J9)) as 

£(/, g) = (-2Cf, g), V / G F 2 ' 2 , V 9 e ff 1 ' 2 

which, compared with JSJ, proves that 2£ to be the unique operator associated 
with the Dirichlet form £ . Further, from self-djointness of £, we infer 

(10) {Lf,g) li = {f,Lg) tl , Vf,g€H 



2.2 



where (-, ■) refers to the usual L 2 inner product. We can take g = 1 to get that \i 
is an invariant measure for the process X t . This proves the claim. □ 
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When is the function exp(— 2k(x)) integrable ? This question is critical to both 
the recurrence property of the diffusion process as well as the existence of a unique 
long term stationary distribution. Its answer, however, is geometric in nature. It is 
intuitive that k needs to be nonnegativc. Notice that if k is nonnegativc, by virtue 
of being positively homogeneous, it is the gauge function (Minkowski functional) of 
a set containing the origin. That is to say, if we define the unit ball and the surface 
given by k respectively as 

(11) C = {x : k{x) < 1}, S = {x : k(x) = 1}, 

it is not difficult to see that k satisfies the relation k{x) = inf {a > : x G aC}. 
We have the following definition. 

Definition 1. A continuous, nonnegative, positively homogeneous function k : 
M. n ->RU {oo} is said to be irreducible if it satisfies k(x) = if and only if x = 0. 

Lemma 5. // k is continuous and irreducible either k(x) > 0, V x G l n , or 
k(x) < 0, V x G R n . Moreover, if k(x) > 0, V a; £ K" , x ^ 0, then Vol (C) < oo. 

Proof. Suppose that there are points xq and x\ such that k(xg) > and k(x\) < 0. 
We can choose a continuous curve 74, t G [0,1], in R n such that 70 = x\ and 
71 = X2 and ^ 7[0, 1]. Since k is continuous, by the intermediate value theorem, 
there exists at* € (0, 1) such that k(j t *) = but 74. 7^ 0. But this is impossible if 
k is irreducible, and we have proved the first assertion of the lemma. 

For the second assertion, we need to show that C is bounded. Suppose, on the 
contrary, we can find a sequence {x n } C C such that limn^oo ||x rl || = 00. One can 
assume that \\x n \\ > 1, for all n € N. Then the points y„ = x n j ||x ra || satisfy 

k(vn) = k(x n )/\\x n \\ < 1, VueN. 

Thus, y n G C, for all n = 1,2,.... However, there exists a subsequence of y n , say 
{Un m } such that linim^oo y nm = z, for some z with ||z|| = 1. Hence, by continuity 
of k, we infer 

k{z) = Urn %„J = lim k j Xnm ] = 0. 

m^oo m^oo \\X n 

The final equality is due to the fact that < k(x n ) < 1 for all n and lim ||a; n || =00. 
Since z ^ 0, this contradicts our assumption that k is irreducible. Hence we are 
done. □ 

We now show that the process in (0 is Harris recurrent if k is nonnegative and 
irreducible. It then follows (see [8], Section 7.5) that it has a unique invariant 
measure /x described above in Proposition [H Moreover, if Pt{x) is the marginal 
distribution of X t when Xq = x, then lim^oo |jP t (x) — A*|| T y = 0. Here ||-|| TV 
refers to the total variation norm on measures. 

The following claim settles the argument. 

Lemma 6. Consider the notations and assumptions in Proposition® Suppose that 
the function k is a nonnegative, irreducible, positively homogeneous function. Then 
the process X t is recurrent. 

Proof. We will use Corollary 7.5.4 in [8]. We need to consider the quantity 

d(x) = n + 2 {x, b(x)} , xeR n . 
By our definition we have (x,b(x)) = —k(x). Thus, d{x) = n — 2k(x). 
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Now, since k is non-negative and irreducible, it is growing to infinity uniformly 
in all directions radially outward from zero. The way to see this is to note 

k(x) = ||ar||A;(x/||x||) > ||x|| inf k(y) = c x ||x|| 

The constant c\ = inf|| y || = i k(y) is positive since k is a strictly positive continuous 
function on the compact set {y : \\y\\ = 1}. 

Now, if we fix an e > 0, there exists R > such that d(x) < — e for all x with 
1 1 a; 1 1 > R. Let Tr be the first hitting time of the compact set Br = {x £ R n : 
\\x\\ < R}. By Corollary 7.5.4 in [8|, we immediately obtain E X (T R ) < \\x\\ 2 /e. 
Thus Br gets visited infinitely often and hence the process X t is recurrent. □ 

The integrability of the function exp(— 2k(x)) requires precisely the same condi- 
tion as in the last lemma. 

Lemma 7. Suppose that k : R" — > R is a nonnegative, irreducible, continuous, 
positively homogeneous function. Then the exp(— 2k(x)) is an integrable function. 

To prove the previous lemma we need the following polar decomposition formula. 
For any nonnegative, irreducible fc, define the surface projection function 8 : R n — > 
S by " 

(12) e(x) = -^ T , x £ M™\{0}, and O(0) = 0. 

k{x) 

Clearly the range of G is the surface S defined in (111H . 
We take the following two measurable spaces 

Mi = (r, B(R n )) and M 2 = (r + x S, B(R+) ® B(S) 

One can construct a measurable map T : Mi — > Mi given by T(x) = (k(x), Q(x)), 
Vie R n . It clearly follows from the definition that T is a one-to-one map. We 
prove the following slightly general result for the sake of completeness. 

Lemma 8. For any nonnegative, irreducible, positively homogeneous function k : 
R" — > R + U {oo}, and any integrable f : R n — » R, we have 



(13) / f{x)dx = nVo\(C) / r"" 1 / f(r-z)dfi(z)dr, 

J {x:k(x)<oo} JO JS 

where dx refers to the n-dimensional Lebesgue measure and fj, is the cone measure 
on S defined as 

(14) ^ ) = _l_Vol(T- 1 ([0,l]x J B)), VEeB(S). 

Proof. We first prove (fl"3|) for functions / equal to indicators of sets A = T -1 ([0, b] x E), 
where E £ B{S) and b > 0. We have 

Vol (A) = J 1a&x = J l{k{x) <b}l{Q(x) £ E}dx. 

If we make let y = b _1 x, then, by positive homogeneity of k, one can write the last 
equation as 

Vol (A) = b n J 1{%) < l}l{9(y) G = 6™Vol (T' 1 ^, 1] x £?)) 

= b n Vo\ (C) ■ n(E) = Vol (C) n [ r™" 1 / d/i dr, 
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which proves (|13|) for this particular case. The rest of the argument follows from 
standard measure theoretic approximation results. □ 

Proof of Lemma^\ Since 2k is another nonncgativc. irreducible, positively homo- 
geneous function, it suffices to show that exp(— k(x)) is intcgrablc. 

By Lemma [5j the set C = {x £ R n : k(x) < 1} has a finite volume, and hence 
the cone measure on S = dC is well-defined. From the change of variable formula 
in Lemma [8] we then obtain 

/>oo 

e -H*)dx = nYol(C) / r"- 1 e^ r dr 



o 



oo 



nVol(C) / s n - 1 e- s ds = Vol (C)nT{n) 
Jo 

Vol(C)r(n+ 1) < oo. 



□ 



Note that, under the probability measure with the unnormalized density e~ klyX \ 
the random variable k(X) is always distributed as Gamma(n) irrespective of k and 
independently of the vector Q(X). Similarly, under the uniform measure on C, the 
random variable k(X) is always distributed as Beta(n, 1) independently of Q(X). 
Under both these measures, Q(X) has the same law. This provides a link between 
the two probability measures which is important in their understanding. 

What can be recovered when k is not irreducible ? In general little, except when 
k is the gauge function of a lower dimensional polytope. The following proposition 
generalizes Proposition 2J 

Proposition 9. Let k : W 1 — > R be a continuous function whose derivative in the 
sense of distribution is represented by a bounded function. Suppose there exists a 
subspace H C R" such that if y £ H and z £ H 1 ^ , then 

k{y + z) = k 1 {y) + k 2 (z). 

Additionally, assume that k\ is a nonnegative, irreducible, positive homogeneous 
function on H. Consider the solution to |[7J) when b(x) represents — Vfc(a;), where 
V is in the sense of distributions. Assume that a solution to with this drift 
exists. Let A : K n — > H be a projection matrix onto the subspace H . Then 

(1) the process Y t — AX t has a unique reversible stationary probability distri- 
bution \x. 

(2) Suppose k' is any other function defined as 

k'(x) = k 1 (Ax) + k' 2 (x-Ax), Vief", 

for some non-negative function k' 2 such that exp(— 2k') is integrable. Then 
/i is the law of the random vector Y — AX , where X is a random vector 
with density proportional to exp(— 2k'). 

Proof. Define k' 2 : R" -> R to be 

(15) k' 2 (x) = sup I, 

l<j<n-d 

where i = 1, . . . ,n — d is an orthogonal basis of H 1 -. Modify the function k by 
defining 

k'(x) = ki(y) + k' 2 (z), V y £ H, z £ H^, x = y + z. 
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• We claim that k' is a non-negative, irreducible, positively homogeneous function 
on M™. Let us verify the condition in definition [T] Since k\ and k' 2 are both positive, 
for x = y + z, y € H, z£ H , we get 

{h'{x) = 0} {k'{x) = 0} {fci(y) = 0} and {^(z) = 0}. 

But, by assumption for fci and by construction for k' 2 , we get that k'(x) = if and 
only if y = and z = 0, that is, if and only if x = 0. This proves our claim. 
Consider the solutions of the following two stochastic differential equations 

(16) dX t = -Vk(X' t )dt + dW t , 

(17) dX' t = -Vk'(X' t )dt + dW t , 

where the drifts are function representing the derivatives in the sense of distribu- 
tions. That X[ exists (in the weak sense) is clear. 

• We claim that X' has a unique reversible invariant distribution fi' with unnor- 
malizcd density exp(— 2k'(x))dx. From Lemma [3 we deduce that exp(— 2k') is 
integrable. Hence we can suitably normalize and get a probability measure whose 
density is proportional to exp(— 2k'{x)). Now, we can apply Proposition 0] to claim 
that the unique reversible invariant distribution of the process X' exists and is 
given by the unnormalized density exp(— 2k'(x)). 

Now, consider a projection matrix A with range space H . Thus, A 1 - = I — A is 
a projection onto H 1 - . Now, by assumption, the function k splits additively as 

k(x) = kx(Ax) + k 2 {A- L x) 

Taking gradient on both sides, one obtains 

Vk(x) = A'Vki(Ax) + (A ± )'Vk 2 (A ± x). 

Here A' and (A- 1 )' refers to taking adjoints. But A being a projection is self-adjoint 
and satisfies A 2 = A. Thus, it follows that 

(18) AVk{x) = A 2 Vfci(Ac) + A(I - A)Vk 2 (A ± x) = AVk^Ax). 
Exactly in the same way we also have 

(19) AVfc'(x) = AA'\7k x {Ax) + A{A^)'\7k' 2 {A^x) = Ak x {Ax). 

Now, consider the processes Y = AX and Y' = AX' which clearly satisfy the 
following differential equations 

dY t = -AVk{X t )dt + AdW t , and 
^ dYl = -AX7k'(X' t )dt + AdW t . 

Using equations (fT8|) . and (fl9|) we can rewrite the above equations as 
dY t = -AVh(Y) + AdW t , and 
dY( = -AVkx(Y') + AdW t . 

Clearly the laws of the process Y and Y' are identical due to the uniqueness in 
law of the weak solutions of their stochastic differential equation. To get past the 
arbitrariness of the hyperplane H, one can simply observe that if dimcnsion(_ff) = 
d < n, there exists a (d x n) matrix D which is a bijection between H and R d . 
The laws of DY and DY' arc identical by standard theory of SDE. Now one simply 
inverts D onto H to obtain our conclusion. 
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Now, since A is a linear map, the process Y' has a unique invariant distribution 
v induced by the invariant distribution of X' ', and given by 

(21) y{B) = n' '{Ax 6B), VBe B(H), 

where \i' is the unique invariant distribution for the process X' . Thus, by the 
equality in law of the processes Y and Y' it follows that the process Y t — AX t has 
a unique invariant distribution given by v according to the recipe above. 

It is easy to give an intrinsic description of v from (|2"Tj) . If X' is distributed 
according to fj,', then for any B € B(H), we have 

v(B) = fjf (AX' € B) = I cxp (-2ki{Ax) - 2k' 2 {x - Ax)) dx 

JAxGB 

= / cxp (— 2k\(y)) dy I cxp(—2k' 2 (z))dz. 

JyeB JH 1 ^ 

Clearly, v has a density proportional to exp(— 2k\(y)) with respect to the Lebesgue 
measure restricted to the hypcrplanc H. □ 



3. SlMPLICIAL CONES AND EXPONENTIAL DISTRIBUTIONS 

We have seen in the last section that the invariant probability distributions for 
the SDEs described in Q have unnormalizcd densities with respect to the Lebesgue 
measure given by exp(— 2k{x)), where is a nonnegative, irreducible, positively ho- 
mogeneous function. Even with such an explicit description of the density function 
it can be very hard to compute any means, variances, or one-dimensional marginal 
distributions. Our objective in this section is to link them to the combinatorial 
structures of the unit ball generated by the positively homogeneous function 2k, 
particularly when these are star-shaped or convex polytopes. In the special case of 
simplicial polytopes this allows us to furnish a complete description of the invariant 
measure in terms of independent Exponential random variables. 

The main geometric idea is the following. Consider, as before, a drift function 
b which is constant over cones Ci, C2, . . . , C r . We assume that b is the negative 
gradient in the sense of distributions of the function k{x) = — (x, b(x)). Thus, 2k(x) 
is a linear function inside each cone Cj. 

Suppose now that d is simplicial. Simplicial cones are those that can be trans- 
formed to the positive quadrant by applying a non-singular linear transformation. 
That is to say, for each C, there exists n linearly independent vectors in K™, de- 
noted by /3i, /?2, ■ ■ ■ , fln, such that d = {x e E™ : x = J2i a iPi, > V i}. It is 
now not hard to see that in that case the probability measure given by exp(— 2k(x)) 
restricted to the cone Ci must be a linear transformation of independent Exponen- 
tial distributions. The difficulty in the execution of the previous argument is to 
identify from the function b that a simplicial polytope is lurking behind the scenes 
and to compute the necessary linear transformations. In the following two subsec- 
tions we consider two general classes of examples where the argument can be fully 
carried out. One, the regular case, is where each Ci can be mapped to any other 
by a group of orthogonal transformations. This leads us to a connection with finite 
irreducible Coxeter groups. The other, which is not regular, deals with graphs and 
the combinatorics hidden in their structure. 
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3.1. The regular case: groups of orthogonal transformations. Consider a 
finite dimensional Euclidean vector space V . A linear transformation from V to 
itself is called orthogonal if the corresponding matrix is orthogonal. In that case 
the determinant of the transformation is ±1. An important class of orthogonal 
matrices is given by reflections. A reflection along a unit vector u G V corresponds 
to the matrix / — 2uu' (all vectors are columns and u' denotes the transpose of 
it). Geometrically it produces the mirror image of any vector with respect to the 
hyperplane orthogonal to u. A reflection group is a group of matrices, each element 
of which is a reflection matrix. 

Let G be a group of orthogonal matrices. G is called irreducible if there is no non- 
trivial subspace W of V which is stable under the action of G, i.e., p(s)W C W, for 
all s G G. As a recurring example, considering the symmetric group of permutations 
on n elements. It has a natural representation as permutation matrices on M™. This 
is not irreducible since the one-dimensional subspace W\ spanned by the vector of 
all ones remains invariant under the action of the group. However, the action 
restricted to W\ and W2 = W^- is irreducible. 

For a finite group of orthogonal linear transformations G we now define what is 
known as a fundamental region. Please see Chapter 3 in |10j for more details. 

Fundamental Region. A subset F C V is known as a fundamental region for a 
group of orthogonal transformations G if 

(1) F is open, 

(2) F n AF = if A I , where I is the identity matrix and A G G, 

(3) V = Ua^gAF, where B denotes the topological closure of a subset B. 

For the representation of the symmetric group as permutation matrices such a 
region is provided by the cone {1 £ R" : x\ < x 2 < ■ ■ ■ < x n }. In general, funda- 
mental regions are not unique. 

Henceforth we will work with V = l n even though proper subspaces of K" 
provide another rich class of examples. 

The groups we will be interested in, which includes the permutations as a special 
case, are generated purely by reflections. When irreducible, these groups are known 
as Coxeter groups and we review their basic structure below. 

A reflection along a vector r is uniquely characterized by the fact that it keeps 
every vector orthogonal to r unchanged and flips the sign of every multiple of r. A 
Coxeter group is a finite irreducible group of orthogonal transformations generated 
by finitely many reflections. A frequent class of examples are the Dihedral groups 
which are the symmetry groups of regular polygons. Although, these groups contain 
both rotations (think of a unit square being rotated by tt/2) and reflections (the 
square getting reflected along the mid- axis), one can show that these groups can 
be generated purely by the reflection elements (reflecting the square twice along 
different axes amounts to a rotation). 

Closely associated with the Coxeter groups is the idea of root systems. If G 
contains a reflection along r, then both r and — r are known as roots of G. Let A 
be the set of all the roots of G, usually referred to as the root system of the Coxeter 
group. Recall the definition of a fundamental region. We are going to construct a 
fundamental region for G which is a simplicial cone. 

Let G be a Coxeter group acting on E™. Thus, in particular, it is irreducible. 
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Choose any vector u £ 1™ such that (u, r) ^ for any root r of G. Then the 
root system is partitioned into two subsets 

A+ = {r e A : (r,u) > 0}, A" = {r e A : (r,u) < 0}. 

Theorem 4.1.7 in [lOj (and several lemmas preceding it) proves the following result. 

Theorem 10. There is a unique collection of n many vectors Tl u in AJ such 
that every vector in AJ can be written as a linear combination of elements in Tl u 
with non-negative coefficients. This collection, known as the u-base or fundamental 
roots, and denoted by W u , is linearly independent and forms a basis o/R™. 

A u-base provides a fundamental region for a Coxeter group G, by defining 

(22) F u := {x £l":(i,r)>0 V r G II„} . 

That this is indeed a fundamental region has been proved in Theorem 4.2.4 in 
[TO] , That this is a simplicial cone follows since the vectors in il u are linearly 
independent. 

Let us now return to the framework in Proposition 3J We start with a drift 
function b : R™ — > R that is scale invariant, i.e., b(ax) = b(x) for all a > 0. We 
would like to analyze the probability measure given by normalizing cxp(— 2k(x)) 
where k(x) = — (x,b(x)). Suppose the drift function b takes finitely many values 
on R™ and satisfies the property that b(Ax) = Ab(x) for all A £ G for some Coxeter 
group G. Can we claim that there is a finite sequence of non-overlapping simplicial 
cones whose union is the whole space and such that b takes a constant value inside 
each cone ? 

The answer is no in general. However, there is a simple sufficient condition 
which indeed guarantees an affirmative answer to the question. Let, as before, 
k(x) = — (x, b(x)}. If k is non-negative, then clearly, by positive homogeneity, k is 
the Minkowski functional of a star-shaped body containing the origin. The question 
in the previous paragraph is equivalent to asking whether this star-shaped body can 
be triangulated as a disjoint union of simplices each of which contains the origin 
as a extreme point. We are going to show below that the answer to the question is 
yes, if k(x) = nvaxAec (A, Ax) for some vector A. Seen through a geometric lens, 
this is equivalent to the statement that the star-shaped body generated by k is a 
convex polytope which then turns out to be simplicial. 

Lemma 11. Let G be a finite irreducible group of orthogonal matrices on R n . For 
any A € R™, A ^ 0, the function 

(23) k(x) = max (A, At) , i£R", 

is a nonnegative positively homogeneous function on l n , which is irreducible. More- 
over, k is invariant under the action of the group. That is k{x) = k (Ax), for all 
A e G, and for all x g R" . 

Proof. It is trivial to see that k is positively homogeneous. To show that it is 
nonnegative we use the fact (see [6]) that for any non-trivial irreducible group G of 
orthogonal matrices, the sum Xmsg ^ 1S ^ ne zer0 matrix. Since G is irreducible, it 
follows that 

(A, Ax) = /\,Y Ax ) = 0- 

AdG \ A I 
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But this implies that the maximum must be non-negative. Thus k is nonnegative. 

To prove that k must be strictly positive for all nonzero vectors, note that for 
the previous argument, k(x) = for some x would imply that 

(X,Ax}=0, VAeG. 

Define the subspace 

(24) V x = {x G R" : (A, Ax) = V A g G} . 

Next, note that V\ is stable under the action of G. To see this take any B £ G and 
any x £ V\, then clearly Bx £ V\. Thus AV\ C V\ for all A £ G. But, since G is 
irreducible, V\ must be either zero or the entire space. If V\ is the entire subspace, 
then by putting A to be the identity in the definition ([24]) . we get (X,x) = OV x £ V. 
This shows that A must be zero, which we have ruled out in our assumption. The 
invariance of k under the action of G is clear by the homomorphism property. □ 

The last lemma proves that k is the gauge function of a convex polytope con- 
taining the origin. The following lemma shows that the number of extremal faces 
of the polytope is given by the size of the orbit of A. 

Lemma 12. Consider k as in Lemma \TT\ Given any A, there exists x ^ such 
that k(x) = (A, x) . Moreover, k (Ax) = (AX, Ax) for all A£ G. 

Proof. To prove the first assertion, suppose that for all x, we have 

(A, x) < k(x) = max MA, x) . 

AeG 

Then, for any such x, for any B £ G, we also have 

(BX, x) = (A, B'x) < k (B'x) = max (AX, B'x) = max (BAX, x) = k(x). 

But that would imply maxAeG (AX, x) < k(x) which is clearly a contradiction. 
For the second assertion, consider x, X such that k(x) = (x,X). Now 

k (Bx) = max MA, Bx) = max (B' AX, x) = max ((B~ 1 A)X, x) . 

AeG AeG AeG x 1 

The right hand side is maximized when B~ 1 A = I which proves the lemma. □ 

We now show that for Coxeter groups that the polytope generated by k is sim- 
plicial. That is to say, all its extremal facets are simplices. 

Lemma 13. Consider any n- dimensional irreducible group of orthogonal matrices 
G. Let X £ R™ be such that AX ^ X for all X ^= e. In other words, X has no 
non-trivial stabilizer. 

Consider the region C = {x £ W l : k(x) = (A, x)}. Then the interior of C , given 
by Co = {(A, x) > (AX, x) , V A £ G}, provides a fundamental region for the group. 
Additionally, if G is a Coxeter group, C is an n-dimensional closed simplicial cone. 

Proof. Note that C is the region {x £ R" : (A, x) > (AX, x) MA £ G}. We first show 
that C is a n-dimcnsional convex cone. Label the non-identity elements of the group 
G by Ax,A 2 , ...,A N where N + 1 = \G\. 
Consider the matrix Q defined by 

X-AxX 
X — A 2 X 



(25) Q = 



X - A N X 
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where all vectors are row vectors. This matrix Q when applied to vectors of C 
produces nonnegative entries. The dimension of Q is N x n. We first show that the 
rank of Q is n. Note that, trivially the rank cannot be more than n. We show that 
the dimension of the kernel is zero which proves that the rank must be exactly n. 

Let K, denote the kernel, {x : Qx = 0}. Then we claim that K, is invariant under 
the action of the group. This is because, x € K, iff (A, x) = (AX, x) for all A e G. 
But, for any B, A G G, we also have 

(AX, Bx) = (B~ 1 AX, x) = (X, x) , since x G K, 

= (B^X, x) = (B'X, x) = (X, Bx) . 

Thus Bx £ K,. This proves that K, stable under the action of the group. But since 
the representation is irreducible, this implies that K, is either zero or the full space. 
But, it is easy to see that if K, is the full space, then A must be zero. This proves 
that C is ra-dimcnsional. That it is a convex cone is obvious. 

Since the dimension of C is n and the stabilizer of A is the identity, the interior 
of the cone is given by Co- We now show that Co is a fundamental region for G by 
verifying the definition. Co is open by definition. For any A £ G, A ^ I, note that 
ACq is the following set 

{Ax, xeC Q } = {y: (A, A- l y) > (BX, A~ l y) , VB G G) 

= {y:(AX,y)>(BX,y),VBeG}. 

Thus x G Co H ACq implies (X,x) > (AX, x) > (X,x) which is impossible. Thus the 
intersection must be empty. It is also trivial to see that UaegAC = M™. This shows 
that Co is a fundamental region. 

For Coxeter groups we now show that Co is the same region as F\ defined in 
(|221) . Notice first that if A is a reflection along a vector r for some A £ G, then 

(26) (X-AX,x)=2 —2 . 

IMI 

Now, suppose x e Co- Then (A — AA, x) > for all non-identity A <G G, in partic- 
ular, for all A which corresponds to reflections along the roots. Thus, for any root 
r e A^, from the above equality we get that (r, x) > 0. From the definition of F\, 
it is now obvious that igFj. Thus we have shown that Co C F\. 

For the reverse equality, note that if Co is a proper subset of Fa , then for every 
A G G, the set ACq is a proper subset of AF\. But, each AF\ is disjoint and the 
union of the closures of ACq is the entire W . This is clearly impossible. Thus, 
we have shown that equality holds among the two fundamental regions Co and F\ . 
Since F\ is a simplicial cone, so is Co- Thus C is a closed simplicial cone. □ 

The connection between simplicial cones and Exponential distributions is made 
precise in the next lemma. 

Lemma 14. Consider a sequence of simplicial cones C\, C2, . . . , C r which are open, 
disjoint, and the closure of their union is the whole space. Let k be a nonnegative, 
irreducible, positively homogeneous function such that k is linear on each Ci. That 
is, k(x) = (A,-, x) , for all x G C,, for some sequence of vectors Ai, A2, . . . , A r which 
may not be all distinct. 

Let X be a random variable whose density with respect to the Lebesgue measure 
on I™ is proportional to e~ 2k<yX ^ . Let Bi, i = 1, 2, . . . ,r, be any set of invertible ma- 
trices such that for each i, the matrix Bi maps the cone Ct onto the n-dimensional 
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quadrant. Also, let a\(i), . . . , ct n (i) be the coefficients in the unique representation 
X[ = a'(i)Bi. Then, the random vector Y = (Y\, Y2, ■ ■ ■ , Y n ), where 

Y 3 := aj (i){(B i ) j „X), ifXed, 

is a vector of iid Exponential(2) random variables. Here (-B;)j* denotes the jth row 
of the matrix Bi . 

Proof. Since the cones Cj's are simplicial, the existence of the matrices B^s follows 
from the definition. Moreover, it follows from the definition of a(i) that if x G C- L 
and z = BiX, then k(x) ~ (Xi,x) = J^ - aj(i)zj irrespective of i. Let B* be the 
matrix given by 

(B*) j , = a j {i)(B i )^ > i = l,2,...,d. 
Since the transformations are piecewise linear, it follows that 

Y^J^B'XliXeCj) 

3 

has a density proportional to exp{— 2 ^ . yj} over the quadrant {y : y\ > 0, j/2 > 
0, ■ ■ ■ ,Ud > 0}. This immediately identifies itself as the joint density of iid expo- 
nentials with rate two. This proves part (2). □ 

Proposition 15. Let G denote a Coxeter group acting on W™. For X <G M. n , X ^ 0, 
let k(x) = max^gG {AX, x) . Let v be the probability measure with unnormalized 
density exp{— 2k(x)} on W 1 . Then the following statements hold true. 

(1) When the stabilizer of X is trivial, the conic hull of the finite set {A — 
AX, A S G} contains n linearly independent generating vectors {rji, 772, ... , r]n}- 
That is, every other vector in the set can be expressed as a linear combina- 
tion of the generators with nonnegative coefficients. 

(2) Let X denote a random variable with distribution v . Also let on denote the 
unique positive coefficient of r)i in the expansion 

n 

X = y^a»?/i- 

i=X 

For i = 1, 2, . . . , n, define the change of variable 

Y t = a, (Arji, X) , when k(X) = (AX, X) . 

Then, the vector Y — {Y\, Y2, . . . ,Y n ) are iid Exponential(2) random vari- 
ables. 

Proof. To prove the first assertion, we use Lemma [131 Assume first that the stabi- 
lizer of A is trivial. Then, Lemma [T31 tells us that the cone 

(27) C = {x: (X-AX,x) > 0} 

is an ri-dimcnsional simplicial cone. Hence there exists exactly n many linearly 
independent generators among the set {A — AX, A £ G} such that every other 
vector is a linear combination with nonnegative coefficients. The rest of the result 
follows directly from Lemma [141 Notice that the coefficients of A in the expansion 
A = 2™=i a i 7 li are positive by Farkas lemma. This is because, for any non-zero 
y € C, the inner product (A, y) = k(y) > by irreducibility and nonnegativity of k. 
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When A has a nontrivial stabilizer, the cone in (|27|) is a union of several simplicial 
cones. The simplest way to see this is to take a sequence A; which have no nontrivial 
stabilizers and which converges to A. The component cones are then given by the 
limits of the sequence of simplicial cones generated by them. In any case Lcmma [T4l 
still holds, however, the vectors {771, . . . , rj n } have to be determined by the limiting 
procedure. □ 

The proof of Proposition [2] in the Introduction now follows easily. 

3.2. Examples. Let us consider some examples of consequences of Proposition [21 



Example 1: Rank based interactions. Brownian motions with rank based 
interactions have been considered in equation ([1]). Clearly the drift function b is 
constant over finitely many cones determined by the permutation generated by the 
ordered coordinates. Let x[l] < x[2] < ...x[n] denote the coordinates of an n 
dimensional vector arranged in increasing order. It is easy to see that 

n 

k(x) = y]Sjx[i], 
i=i 

is a positively homogeneous function which is not irreducible since it takes a con- 
stant value over the linear span of the vector of all ones. However, if we let H be 
the subspace orthogonal to 1, then k splits additively as 

n 

k{x) = — Si (x[i] — x) + ndx. 

i=l 

Let fci : H — ► ffi denote the restriction of k to H, then it is clear that ki(x) = 
ki{A a x) for any permutation matrix A a . Now, the group of permutation matrices 
acting on H is well known to be irreducible and generated by reflections along 
Ci + \ — e.i for i = 1, 2, . . . , n — 1. This is just a restatement of the fact that every 
permutation can be written as a product of transposes. Thus, it is a Coxeter group, 
often denoted by A n -i- Thus the conclusions of Proposition [TT1 applies and fei is 
irreducible if 

(28) fci(x) = max (-A a S, x) . 

This condition is equivalent to the condition b\ > 82 >•••_: 6 n . 

Hence, from Proposition [9] it also follows that the projection PnXt of the diffu- 
sion on to H has an invariant distribution whose density with respect to Lebcsgue 
measure on H is proportional to exp(— 2ki(x)). 

We now apply Proposition [T5J One can see that the set of vectors —5 + A a S, as 
a ranges over permutations, contains positive multiples of vectors e^+i — e^, since 
they correspond to the transposition of i and i + 1 . These n—1 linearly independent 
vectors are the conic extremes of the set. Thus, by Proposition [T"5l the spacings 
X [i + 1] — X [i] arc independent Exponential random variables under the invariant 
distribution. The correct rates can be easily verified. 

In this example it is easy to see the shortfall of the sufficient condition (|28p . The 
drift function is constant over a fundamental region F = {x G H : x\ < X2 < 
... < x n } which is clearly simplicial. Now k(x) is irreducible if and only if the unit 
ball generated by k is compact. By symmetry, we can restrict our attention to F. 
Since F is simplicial we can apply a suitable linear transformation to map it to 
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the positive quadrant. Thus, it can be easily verified that if S denoted the average 
of the coordinates of 6, the intersection of the unit ball with F is compact if and 
only if 

k 

a* 5) > 0, Vi=l,2,...,n-1. 

i=l 

This is precisely the condition derived by Pal & Pitman in |13j using the theory 
of reflected Brownian motions and is weaker than the sufficient condition that the 
coordinates of 5 decreases with increasing values of the coordinates of x. 

Example 2: Sign-rank based interactions. 

An example of interactions similar to rank-based can be generated by allowing 
both the rank and signs coordinates to determine the drift. As before, we start 
with the n-dimensional SDE: 

(29) dX t = b(X t )dt + dW u 

where Wt is an n-dimensional Brownian motion. Suppose that the drift function 
takes finitely many values, is scale invariant, and b(Ax) = Ab(x) whenever A is 
cither a permutation matrix or a diagonal matrix with each diagonal entry being 
plus or minus one. Thus, not only that the values of the drift get permuted whenever 
the coordinates get permuted, but also the sign of the drift changes with the sign 
of the corresponding coordinate. 

The group generated by the collection of permutation matrices and the diagonal 
matrices of sign flips is a Coxeter group denoted by B n . Please see pages 66-71 of 
[TO] for more details. We can safely apply Proposition O Thus, the n dimensional 
process under such a sign-rank based interaction is recurrent if there is a vector 
A 6 K", A 7^ 0, such that 

k{x) = — (x, b(x)) = ™-a^ (AX, x) . 

If we restrict the above condition to the cone {x : < x\ < X2 < ■ ■ ■ < x n } we see 
that the vector of drifts b(x) = — A where A satisfis that < Ai < A2 < . . . < A„. 

When this does hold true, X has an unique long term stationary distribution. 
To find the decomposition of this probability distribution in terms of independent 
Exponentials we consider a S all of whose coordinates are non-zero and distinct. 
That is, it has a trivial stabilizer subgroup in B n . Consider the conic hull generated 
by the set of vectors {A — AX, A G £>„}. As in the case of rank-based interactions 
one can see that the generators of the conic hull are positive multiples of the vectors 
e\ and {e^+i — e», i = 1, 2, . . . n — 1}. 

Now we apply the final conclusion of Prop [2j To get the vector of Exponentials 
under the invariant distribution, note that 

n 

fc(*) = 5>MH, 
j=i 

where \x\ [1] < \x\ [2] < ... |x| [n] are the ordered values of the absolute values of 
the coordinates (\xi\ , [xs] \x n \). 

Thus, from Proposition [2] it follows that the random vector (\X \ [1], \X\ [i + 1] — 
\X\ [i], i = 1, 2, . . . , n — 1) arc distributed as independent Exponentials. 

To compute the rates of these Exponentially distributed random variables, one 
needs to compute the coefficient of A with respect to the basis vector e\ and {ei+i — 
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ei, i = 1,2, .. .ri — 1}. A simple computation leads us to the conclusion that the 
corresponding vector of rates of these Exponentials are given by 

' n 

2^A S , 5 = 1,2,..., 



Example 3: Constrained sign-rank based interactions. 

This is a interesting class of constrained sign-rank based interactions where not 
all sign changes of coordinates affect the drift vector. Consider again the stochastic 
differential equation (|29| . Suppose that the drift function b takes finitely many 
values, is scale invariant, and b(Ax) — Ab(x) for all permutation matrices A and 
all diagonal matrices with each diagonal entry being positive or negative one with 
the constraint that there are even number of negative ones. Thus, the sign of the 
drift vector changes when either the ordering of coordinates change or when pairs 
of coordinates have flipped their signs. 

The groups generated by the permutation matrices and the diagonal matrices 
with even number of sign flips is clearly a subgroup of the B n . They form, in fact, 
a family of Coxeter groups, usually denoted by T> n where n denotes the dimension 
of the underlying space. We again refer the reader to pages 66-71 of [10] for more 
details. 

We apply Proposition [2] Thus, the n dimensional process under such a con- 
strained sign-rank based interaction is recurrent if there is a vector A £ IR ra , A =/= 0, 
such that 

k(x) = — (x, b(x)) = max (AX, x) . 

The above condition is more difficult to analyze than the previous examples. One 
can show using known results about the fundamental root systems of V n (page 71 
in [10]) that the drift is determined by the fact that over the cone 

(30) {x : < x\ + X2, xi < x 2 < ■ ■ ■ < x n } , 

the drift is a constant b(x) = —A where A satisfis that < Ai + A 2 and Ai < A2 < 
• • ■ < A ra . This cone is actually a fundamental region for the group T> n . Thus, the 
drift vector is now determined over entire 1™ by the property b(Ax) = Ab(x) for 
all A G T> n . 

Under this condition the process has a long term stationary distribution. To find 
what functions turn out to be independent Exponentials, we need to understand, 
for a given x, what unique A S T> n will produce k(x) = (X,Ax). Clearly, this will 
happen if A is chosen such that Ax belongs to the cone ([3T)]) . There are two cases 
to consider. One, when the number of coordinates of x that are negative is even. 
In this case, one simply flips the signs of these coordinates, and then ranks the 
absolute values to get a vector in ([30]) . Both these actions are permissible since 
they correspond to multiplication by matrices in T> n . The other case is when x has 
odd number of negative coordinates. First, one has to flip the sign of all the negative 
coordinates except the least negative one and then rank all the coordinates. In this 
ordering, cither the absolute value of the second least negative coordinate is less 
than the least positive coordinate in which case, the resulting vector is in ([30| . Or, 
it is not, in which case we need to compare the least negative coordinate with the 
least positive coordinate. Their sum is either positive or negative, and we make the 
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appropriate (zero or two) sign flips to get the right transformation. Let H{x) be 
the resulting vector produced by the above procedure. 

Under this stationary distribution, the vector of random vector H(X\ , X2 , . . . , X n ) 
is distributed as n independent Exponentials. Furthermore, as in the case of sign- 
rank interactions, one can work out the linear algebra to compute that the corre- 
sponding vector of rates of these Exponentials are given by 

(n n n 

$^A», -Ax + ^A S , 2^A S , j = 3,4,...,n 
s— 1 s— 2 s—j 

3.3. An example of irregular interaction. In this subsection we consider an 
example of interacting Brownian motions with drift functions that are still piecewise 
constant on cones, but are not consistent with any group action. 

Consider a graph Q on n vertices where the vetices are labeled by {1, 2, . . . , n}. 
The edge between i and j have an associated edge weight /3ij , which is zero if there 
is no edge between the two vertices. 

Consider the SDE on K n given by 

n 

(31) dX t (i) = £ /%sign (X t (j) - X t (i)) dt + dW t (i), 

i=i 

where, as before, W = (M^(l), W(2), . . . , W(n)) is an n-dimcnsional Brownian mo- 
tion. 

When all the edge weights are nonnegativc, the model can be described by saying 
that the Brownian motions, which are indexed by the vertices of the graphs, get 
attracted towards one another. The constants pij measure the strength of their 
attraction. 

The appropriately defined drift function b(x) is piecewise constant on the family 
of cones CV := {x : x w n) < x n (2) < ■ • ■ < x n r n )}, where tt ranges over all permuta- 
tions of n labels. However, it might not satisfy the condition that b(A v x) = A^b^x) 
where A v is the permutation matrix corresponding to tt. 

Note that the drift function b is the negative of the gradient of the positively 
homogeneous function 

k(x) = ^T fa \xi-Xj\, Vie R n . 

i<j 

It can be easily verified (and intuitive) that if Pij's are nonnegative and Q is con- 
nected, the function k is irreducible when restricted to the subspace H orthogonal to 
the vector of all ones (which we denote by 1). Then the conclusions of Proposition^ 
applies. In particular, if we define 

k'(x) = k(x) + \{x,l)\ 

then k' is integrable and both the probability measures induced by functions exp(— 2k) 
and exp(— 2k') on H must be the same. 

Assume > and Q is connected. For convenience absorb the factor of two in 
exp(— 2k) in the definition of (3. Let [i n be the probability measure whose unnor- 
malizcd density is given by 

(32) exp - 2J Pij \xi —Xj\ — 
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What properties of the probability measure can we explicitly describe ? Clearly, 
any deep inspection of such a general family is extremely difficult. We will improve 
our chances if wc restrict the edge weights to the following class. Consider n positive 
constants to i, m2, Let (3ij = mimj for all pairs i, j. In particular the graph 
is complete. One can think of rrii as the mass of the ith particle, and hence the 
strength of attraction between particles i and j is proportional to the product of 
their masses. In fact, the case when all the to^'s are equal to one have been dealt 
with in [3] where they were named the one-dimensional gravity model. 
Now suppose X = {X\, . . . , X n ) follows the p.d.f. 



C„ exp - 2J rmmj \xi - Xj\ 



i<3 



J2- 



where C„ is the normalizing constant. Let M = Yl7=i mi ^ e ^ ne total mass of 
the system. For each i, let Y, = Xus and il(i) be the (random) index j such that 
Yi = Xj. The joint p.m.f. of (II, Y) at a point (ir,y) is 



C n exp - ^ m^ i m^ j (yj - y t ) 



E 



where 7r = {tt\, . . . , 7r„) is a permutation of {1, . . . , n} and yi < ■ ■ ■ < y n S M. Now 
let Ai = Yi + i —Yi,i = l,...,n — 1. For each i, and each n £ S n , let 



M 



Then M- 1 E]= l+ i ™* 
at a point (7r, S, yi) is 



I-^(tt), andF n (7r) = 1. The joint density of (n, A, Yi) 



exp - ^2 m -x t m -K j \ ^2&k \ - 



h—i 



i=2 j=l 



C n exp -A/ 2 J2 f<W(l - FifrfiSi 



i=i j=i 



nyi 



Now we can easily integrate out Y\ to get the joint density of (II, A) 



Cexp(-M 2 ^F l (7r)(l - i^))<5, 



where C is now a different constant. Thus, conditional on II = w, Yi, . . . , Y n are 
independent, with Yj ~ Exp(Af 2 Fi(ir)(l — Fi(ir))). It is easy to see from this 
observation that the marginal p.m.f. of II must be 

n-l 



C(m)\ 



I 



'lF i ( 7 r)(l- J F i ( 7 r)) : 

where tt = (vri, . . . , n n ) is any permutation of {1, . . . , n} and C{m) is the normaliz- 
ing constant that depends on the values of mi, ... , to„. 

If mi = m-2 = ■ ■ ■ = to„, then this is the uniform distribution on S n . Otherwise, 
it is a non-uniform distribution on the set of permutations. Thus the cost we 
pay for losing the regularity of group actions is that the spacings between the 
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order statistics are only conditionally independent Exponentials, conditioned on 
the random permutation generated by the ranks. 

It is very difficult to see what sort of distributions on the space of permutations 
the probability mass function (|33[) induces. Clearly the p.m.f. is large when Fi(jr) 
is close to zero or one for most values of i. The intuition from gravity predicts that 
heavier particles should be close and should avoid being too high or too low in rank. 

We now show this to be true in a particularly simple case when there is a single 
distinguished particle. Suppose that mi = a and m 2 = m 3 = . . . = m n = 1. We 
consider the joint distribution as before 



^ = C„ cxp ( - ^2 rriimj \xj - Xj\ + |y\ 



i<3 



We are interested in the derived joint distribution of the ranks of each particle given 
by the general expression in the previous section. 

Let us compute the distribution of the rank of the first particle which is distin- 
guished from the others due to a different mass. 

P(rank of the first particle = j) = P(H(j) = 1) 
(34, _ cfi ' . 



11 fi(<7)(l 



where C is a constant depending on n and a. In the following text, we will freely 
use C as the normalizing constant keeping in mind that the constants might be 
different from one another although they only depend on n and a. 

Now, there are (n — 1)! many permutations a such that <r(j) — 1. For any of 
them 



Fi(a) 



a+n — 1 J ' 

otherwise. 



a+i-l 
a+n— 1 



And thus we can rewrite ([M]) as 
(35) 

p(n(j) = i) = c 



nj-l i ( i _ i \ TT"- 1 a+i-1 / 1 _ a+i-l \ 
i=l Q+n-1 I o+n-1 J lli=j a+n— 1 ^ a+n-1 ^ 

C- 



1 



c? - i)!(n - jV- ns(" + n - i - 1) n,r/(a + * - 1) 

n — l\ a(a + 1) . . . (a + n — j — 1) a(a + 1) . . . (a + j — 2) 



C , 

J -IJ a(a + 1) . . . (a + n - 2) a(a + 1) . . . (a + n - 2) 
= C (J ^ J J a(a + 1) . . . (a + n - j - 1) a(a + 1) . . . (a + j - 2). 

We can immediately infer from the previous expression the following fact: 

P(rank of the first particle = j + 1) _ P(U(j + 1) = 1) 
P(rank of the first particle = j) ~ P(n(j) = 1) 

= {n-j)(a + j - 1) 
j(a + n-j - 1) 
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In other words 

P(rank of the first particle = j + 1) > P(rank of the first particle = j) 

iff (n-j)(a + j- 1) > j(a + n - j - 1), 

iff n(a - I) + (n - a + f )j - j 2 > (a + n - 1 )j - j 2 , 

iff 2(a-l)j < n(a-l). 

Thus, if a is more than 1, the probability of the rank being j increases till j = \n/2\ , 
and then strictly decreases. Clearly, the most likely position for the heavier particle 
is going to be the median. On the other hand, if a < 1 , just the opposite happens, 
and we are likely to see the lighter particle either at the top or trailing behind. 

The probability computed in (|35|) . although seemingly unfriendly, is a very fa- 
miliar object. Consider a Polya's urn scheme which has a red balls and a black 
balls. We play a game where at each step we pick a ball at random and replace 
it in the urn with a ball of the same color. It is well known (see Feller [9]) that if 
we play this game for n — 1 steps the probability we pick exactly j — 1 red balls is 
given by 



n 



- l\ a(a + 1) . . . (a + n - j - 1) a(a + 1 ) . . . (a + j - 2) 



j-ly 2a(2a+l)(2a + 2)...(2a + n-2) 



If we compare the previous expression with (|35|) . the differences are merely in the 
expression of the normalizing constants. Thus, if er(l) = II _1 (f) is the rank of the 
first particle, it is clear that a(l ) — I has the same distribution as the number of 
red balls picked in a Polya's urn scheme run for n — 1 steps. 



Proposition 16. For any a > 0, consider the SDE pip with a distribution of mass 
such that the mass of the first particle being a and the rest being of mass 1 . Then 
the sequence of random variables o~(l)/n converges weakly to the Beta(2a,2a). 

Proof. The proof follows from known results about Polya's urn. The factor of 2 is 
due to the fact that we had earlier absorbed the 2 in (1521). □ 
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